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Abstract We show that the linear viscoelastic materials, and more generally the physical 
phenomena to which Biot's relaxation theory is relevant, can be put in correspondance with 
the laws of processes with independent increments. 

In the one dimensional case this correspondence is one to one with subordinators and gives 
rise naturally to a conjugation relation on subordinators. 
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PRESENTATION 

Certaines equations de la physique et les problemes aux limites associes admettent une 
interpretation probabiliste. C'est le cas de l'equation de la chaleur associee au mouvement 
brownien, ou des equations elliptiques lineaires du second ordre associees aux processus de dif- 
fusion, ou encore des operateurs lineaires integro-differentiels verifiant le principe du maximum 
positif auxquels sont associes des processus de Markov avec sauts. Ceci est etudie par la theorie 
probabiliste du potentiel (cf par exemple [4], [11], [12]). Le cas discontinu est une extension 
mathematique du cas continu et correspond a des problemes non locaux plus rarement ren- 
contres par l'ingenieur. De nombreuses autres equations (cf [14]) possedent aussi maintenant 
des interpretations probabilistes. 

Nous montrons ici que les phenomenes viscoelastiques lineaires sont susceptibles d'une in- 
terpretation probabiliste par des processus discontinus. Cette interpretation est analogue aux 
precedentes, en ce sens, que les grandeurs physiques y apparaissent comme l'esperance de 
fonctionnelles du processus. Neanmoins, ici le temps qui regit le mouvement de la particule 
probabiliste n'est pas le temps qui regit le phenomene physique (comme c'est le cas lorsqu'on 
associe au semigroupe de la chaleur un mouvement brownien). Par ailleurs, cette correspon- 
dance peut etre faite de deux fagons qui sont duales et ceci conduit a degager la notion de 
couple conjugue. La presente redaction est une version detaillee de [6]. 

Une fonction <fi : 1R + — > IR + est une fonction de Bernstein si <p est C°°, <p > et (— l) n D n <j) < 
Vn > 1. Ces fonctions interviennent dans la caracterisation des semi-groupes de convolution 
de probabilites sur IR + avec lesquels elles sont en correspondance biunivoque (cf [1]) ainsi done 
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qu'avec les processus a accroissements independants croissants ou subordinateurs. Ces fonctions 
forment un cone stable par composition dont la structure est assez complexe (cf [7], [8]). Elles 
interviennent aussi, ceci est correlatif, dans la theorie du calcul symbolique sur les generateurs 
infinitesimaux de semi-groupes de Markov (cf [1], [15], [18]). 

II est remarquable que ce soit exactement cette classe de fonctions qui intervienne dans les 
phenomenes de relaxation tels que la viscoelasticite lineaire qui sont regis par une mathematique 
a priori sans rapport avec les objets precedents mais qui se trouve mise en connexion avec eux 
par le raisonnement par variables cachees. Plus precisement: Les phenomenes de relaxation tels 
que la viscoelasticite lineaire sont decrits par des operateurs dissipatifs et sont done lies d'apres le 
theoreme de Lumer-Phillips ([23] p. 250) a des semi-groupes a contraction mais ces semi-groupes 
n'ont aucune interpretation probabiliste en general car ils n'operent pas positivement sur les 
fonctions; dans cette situation e'est le raisonnement par variables cachees qui fait apparaitre 
des objets lies a des semi-groupes de Markov. 

1 Thermo dynamique et viscoelasticite 
1.1 Notations 

Les modeles visco-elastiques les plus simples sont les modeles de Maxwell (pour un liquide 
visco-elastique) : 



Lorsqu'un tel systeme est soumis a une force dependant du temps Q(t), la deformation q(t) 
qui en resulte depend lineairement de Q et si les proprietes du corps ne varient pas avec le 
temps cette transformation est connue si on connait la reponse pour 1' impulsion unite appelee 
reponse impulsionnelle f(t) du systeme dont la reponse generale est alors 




et de Kelvin- Voigt (pour un solide visco-elastique): 





Pour un modele de Maxwell on a 
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ou G est la raideur du ressort et r] la viscosite de 1' amort isseur, et done 

/(*) = (^ + ^)W 

Par ailleurs on peut au contraire imposer une deformation au systeme et mesurer la force qui 
assure l'equilibre. Pour une deformation qui passe de a 1 a l'instant puis ne varie plus, la 
force mesuree est la fonction de relaxation r(t) du systeme dont la reponse generate est alors 

Q(t) = j [Q ^ r(t - r) dq(r) 

-—t 

pour le modele de Maxwell r(t) — Ge *> l{t>o}- 

1.2 Limitations imposees par la thermodynamique 

Les restrictions que la thermodynamique des phenomenes irreversibles impose aux phenomenes 
viscoelastiques ont ete etudies et degages par plusieurs auteurs parfois independamment (cf 
[13] [21] [2] [9] [16] [20]). Nous nous plagons sous les hypotheses de la theorie de Biot [2] et 
de meme que celle-ci les considerations qui suivent peuvent s'appliquer a d'autres phenomenes 
physiques (electriques ou chimiques) pourvu que les hypotheses (Principe d'Onsager, existence 
de variables normales, linearite) soient acceptables au moins en premiere approximation. Pour 
la commodite, nous emploierons le langage de la mecanique des solides (cf [20]). 

Soit un systeme soumis a des forces generalisees Qi, i — 1, • • • , n et decrit par des parametres 
geometriques associes qi, i — 1, • • • , n tels que le travail des forces exterieures, s'ecrive J27=i Qt^Qi- 

Au voisinage d'une position d'equilibre stable ou Ton a pris = Vi, le potentiel thermo- 
dynamique du systeme s'ecrit : 

w = \ Yl a am 

Z id 

ou la matrice (a^) est symetrique semi-definie positive (le terme stable est pris au sens large). 
En calculant la variation d'entropie durant un court intervalle de temps (cf [16] chapitre 13) 
on obtient, pour des vitesses q^ supposees petites et en negligeant les forces d'inertie, que la 
puissance dissipee peut s'ecrire : 

ou la matrice (6^) est symetrique d'apres le principe d'Onsager et semi-definie positive par le 
second principe de la thermodynamique. L'equation devolution est alors: 

3D 8W „ 

dq { dqi 

soit 

(1) O, Y."U'I.i ■ h u ( h 

j 

Le cas de la viscoelasticite lineaire est celui ou les coefficients (a^) et (fr^) sont constants. 
Dans ce cas, la relation lineaire entre les histoires des forces (Qi(t)) et celles des parametres 
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(qj(t)) commute avec les translations du temps et est connue par la reponse impulsionnelle 
fijit) du parametre % a la force Qj. 
Ecrivons l'equation ((TJ) sous la forme 

(2) Aq + Bq = Q 

ou A et B sont des matrices n x n symetriques semi-defmies positives et supposons d'abord 
que B est definie positive. Considerons sur H n la structure euclidienne associee a B dont le 
produit scalaire est 

(u, v) b =< u, Bv >= l uBv 

(<.,.> etant le produit scalaire usuel sur IR n ). Pour cette nouvelle structure euclidienne 
l'operateur B~ X A est auto-adjoint: 

(u, B~ x Av)b =< u,Av >=< Au,v >= (v,B~ 1 Au) B 
Done il existe une base £>-orthonormale (ipk)k=i,...,n sur laquelle B~ x A est diagonale, e'est a dire 

{^j,B- l A^ k )B = <^ j} Aifj k > = 
{^j^k)B = <ipj,Bip k > = 

en particulier 

Aifj k - \ k Bip k = Vfc. 

Soit q(9) [resp. Q(9)] la transformee de Laplace de q(t) [resp. Q(t)] (q(9) = e~ 6t q{t)dt). 
L'equation (J2J) donne 

(3) (A + 9B)q = Q. 
Done, si q{6) est representee sur la base (i^ k ), 

n 

(4) = E&W* 

fc=i 

on a 

(X k + 9)£ k = (Q,il> k ,) B . 

Si nous supposons maintenant que seulement m des n parametres (m < n) sont observables, 
i.e. if Q = (Qi, • • • , Qm, 0, . . . , 0) il vient 



ilj^k 
A fc >0 if ^ = Ac 



if j ^ fc 

1 ifj = k 
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et par (4) 



m n 1 

a = 12 QA12 \ , fl Wj^l- 
,■=1 fc=i A fe + y 

Ainsi, designant par fij(t) the reponse de a l'impulsion unite de Qj, on a 

n 1 

•4 = 12 \ .n ^kjW^jlflki 
fc=l A fc + ^ 

et done 



k=i J0 



On voit facilement que si des hypotheses plus larges sont prises sur la matrice B autorisant 
des valeurs propres nulles ou infinies, on obtient la forme plus generate 

n 

(5) / y (t) = £(1 - e- Afct )4 } + tiy + 

fc=i 

ou les matrices (J^)k = 1, . . . , n, L, K sont symetriques semi-defmies positives. 

Partant de l'equation (5) et passant a la limite on obtient que les reponses impulsionnelles 
des systemes visco-elastiques sont 
dans le cas d'un seul parametre observable 

(6) fit) = ( (1 - e- xt )v(d\) + bt + c 

ou v est une mesure a-finie positive on JR* + =]0,oo[ telle que J °° Y^dv{x) < +oo, et b > 0, 
c > 0; 

et dans le cas de m parametres observables de la forme 

(7) fij(t) = I (1 - e" At )^-(dA) + tLn + Kn, 

ou v = (i/ij) est une matrice semi-definie positive de mesures cr-finies sur satisfaisant 

r x 

/ — — d\uij\{x) < +oo = l,...,n 

et ou les matrices L and K sont symetriques semi-definies positives. 

Ce passage a la limite peut se faire selon des arguments physiques, en prenant l'ensemble 
des limites simples des fonctions de la forme (5), ou en considerant un continu visco-elastique, 
les matrices A et B etant remplacees par des operateurs auto-adjoints non bornes sur un espace 
de Hilbert. Alors l'argument ci-dessus s'etend et donne directement les formes (6) et (7) par la 
representation spectrale des operateurs auto-adjoints. 

Reciproquement, les systemes caracterises par (7) verifient les exigences de la thermody- 
namique et de la stabilite (cf [20] annexe XXI-20 ou [16] chapitre 13). 
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Remarque. Dans le cas d'un seul parametre observable (6) montre que la reponse impulsion- 
nelle generate est une fonction de Bernstein. II est done faux de croire ( comme il est ecrit dans 
certains manuels) que le groupement en parallele d'une suite infinie de modeles de Maxwell 
(ou le groupement en serie d'une suite infinie de modeles de Kelvin-Voigt) donne le modele 
rheologique visco-elastique general, de meme que les mesures discretes ne sont pas les mesures 
generales. On perd en particulier tous les modeles ou la mesure v est absolument continue. 
Ainsi, de nombreux modeles sont analytiquement calculables qui ne correspondent pas a des 
groupements d'amortisseurs et de ressorts. 



2 Processus a Accroissements Independants Stationnaires 

Pour la commodite de la suite, tous les processus aleatoires seront indexes par des lettres 
grecques. 

Nous dirons qu'un processus (Y T ) T > defini sur un espace de probabilite (Q,A, P) est un 
P.A.I.S. s'il est continu a droite et tel que Wi < r 2 < • ■ ■ < r n les variables Y , Y T1 — Y , ■ ■ ■ , Y Tn — 
Y Tn _ 1 sont independantes et si la loi de Y T — Y a ne depend que de r — a. 

On note T T = a{Y a ,a < r). Un processus (Y T ) defini sur (Q,A, P), continu a droite, a 
valeurs P, est un P.A.I.S. si et seulement si (formule de Levy-Khintchine) (cf [15] [19]). 

Vu G P m ,Vr > a > 

(8) JE[ex P {i<u,Y T -Y a >}\F a ] 

= exp{(r - < u, > +i < u, n > + J ( e l<u ' x> - 1 - i < u, x > l{\ x \< 1 })du(x)]} 

ou S est une matrice m x m symetrique semi-defmie positive, n e P m , et v est une mesure 
positive sur P m verifiant ^{0} = et / 1 A \x\ 2 dv(x) < +oo. {y est done cr-finie sur P m \{0}). 

Dans cette formule, le triplet (S,??, v) est determine de fagon unique et s'appelle les car- 
acteristiques locales de Y . Notons, cependant, que d'autres formes peuvent etre donnees a la 
formule de Levy-Khintchine (cf [15]) et que le coefficient n s'en trouve modifie. En revanche, 
£ et v ont une signification intrinseque : Y est une semi-martingale pour (JF CT ) et sa partie 
martingale continue Y c (nulle en zero) a pour crochet matriciel : 

(9) <r c ,y c > r =rs, 

et la mesure de Levy v est determinee par les sauts de Y : la mesure dr x dv sur P + x P m est 
la projection previsible duale de la mesure aleatoire : 

fi{dr, dx) = ^{Y^Y a -}5(a A Y a ){dr, dx). 

(T>0 

ou AY a = Y a — Y a _. C'est-a-dire qu'on a la formule: 

(10) e E H ("> °, AY °) = E / dr / H (", r > v)Mv) 

(T>0 J0 JUrn 

pour tout processus H previsible positif. 
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Un cas particulier important est celui des P.A.I.S. a valeurs reelles positives, ou subordina- 
teurs. lis sont caracterises par (cf par exemple [5]). 

VA > 0, r > o > 
(11) Me- X ^- Y ^\f a ] = e - {T - aMX) 



<p(\) = b\ + J {I- e~ Xx )dv{x) 

b E IR+, v mesure positive sur IR + telle que /(l A x)dv{x) < +oo. 

Les fonctions ip de la forme indiquee en ()1 lj) sont les fonctions de Bernstein nulles en zero, 
elles sont en correspondance bi-univoque avec les semi-groupes de convolution de probabilites 
sur ]R + , continus pour la convergence etroite (cf [1]). 



3 Viscoelasticite et P.A.I.S. 
3.1 Cas d'un parametre observable 

Un materiau viscoelastique lineaire sans vieillissement ne presentant qu'un parametre observ- 
able q a pour reponse impulsionnelle (formule (6)) : 

(12) f(t) = L + Kt + [(1- e- tx )du{x) 

ou L et K sont positifs et v est une mesure de Levy de subordinateur. 

On definit une bijection si on associe a ce materiau le subordinateur X tel que Xo = L et 
dont la fonction de Bernstein est f(t) — L (on pourrait aussi considerer un subordinateur a 
duree de vie exponentielle pour traiter la constante L). 

Remarquons que si on se donne un P.A.I.S. Y reel quelconque, Y est une semi-martingale 
pour sa filtration naturelle. Soit Y c sa partie martingale continue. On peut associer a Y un 
materiau visco-elastique en posant Yq = L et 

(13) /(f) = Y 2 + JE[tYf + £ (1 - e- tAY °)} 

0<a<l 



3.2 Cas de plusieurs parametres observables 

Considerons un P.A.I.S. Y a valeur H m tel que Yq soit une constante de H m . On note (£, r], m) 
les caracteristiques de Y (formule (8)). Y est une semi- martingale pour sa filtration naturelle 
, soit Y c sa partie martingale continue. 

On peut associer alors a Y le materiau viscoelastique a m parametres observables dont les 
reponses impulsionnelles sont 

(14) f^t) = Yfri + K[tY«Y« + £ (1 - e^ AY ^ ' 



0<cr<l 

En effet cela donne 
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(15) Ut) = Yfri + fE y + / (1 - e-^l 2 )^ dm(y) 

JR \{o} \y\ 

et les matrices (i/y) ou est l'image par y — > |y| 2 de la mesure • m(dy) sont toutes les 
matrices symetriques semi-defmies positives de mesures sur verifiant / R » (x A l)d\uij\(x) < 
+oo. 

Remarque Notons que, si le materiau 1 [resp. 2] est associe au processus [resp. Y^] et 
si 

y(i) et y( 2 > sont pris independants, le materiau dont le modele rheologique est la mise en 
serie des modeles des materiaux 1 et 2 (c'est-a-dire qui a les memes parametres observables et 
tel que 

fij{t) = fh(t) + fij(t)) es ^ associe au processus Y^ 1 ' + Y^ 2 \ 



4 Materiaux conjugues et dictionnaires 
4.1 Etude des fonctions de relaxation 

Pour voir a quoi correspond la mise en parallele des modeles rheologiques, considerons les 
fonctions de relaxation, en nous plagant d'abord dans le cas d'un seul parametre observable : 

Considerees comme distributions sur 1R nulles pour t < 0, la reponse impulsionnelle / et la 
fonction de relaxation sont liees par la relation : 

(16) f*r' = 5 . 

ou * designe le produit de convolution. Notant " la transformation de Laplace, la formule (6) 
donne : 

(17) (F)(0) = L + ^+[ -^du(x) 6>0. 
Ceci nous conduit a la definition suivante : 

Definition 1 Une fonction h de IR+ dans JR est appelee transformee de Stieltjes s'il existe 
a > et une mesure \i positive sur 1R + tels que : 



W6 > h{9) = a + J -^—d^{x 



e + x 

Si h est une transformee de Stieltjes le couple (a, fi) est uniquement determine (cf [18]). La 
relation (|T7|) exprime que (/'Test transformee de Stieltjes du couple (L, K5 + x ■ u(x)). 

D'apres le theoreme 2 de [18], si h est une transformee de Stieltjes non nulle, il en est de 
meme de jrnm- H existe done un couple (a, fj) unique non nul tel que : 

1 f 1 

a + dfi(x) 



(f')Ce) ^ + 9 + x 



avec a > et / R+ j^d(j,(x) < +oo. 
Par ([To]) , on a done : 
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(r')(0) = ^=J— = a + / ——dp{x) 
(f)(9) Jul + 9x rv y 



d'ou Fori tire : 



/ POO g x \ 

r(t) = a + //{0}5 + / dp(x) 1 

V J10 X / 



{t>0}- 



Ce qui etablit la : 

Proposition 2 Les fonctions de relaxation sont de la forme 

r(t) = a + (38 + [ e~ tx dp(x) 



oil a, (3 > et p est ime mesure positive sur IR^_ teWe gtze J" R * y^:^p(x) < +oo. 

En effet, il suffit de poser p(dy) = yp(dy) ou p est l'image de p\u* + par x — > -. La condition 
/ j^dp(x) < +oo impose alors : 

(18) / -^—dp(y) < +oo. 



1 + y 

Remarque. Ainsi les restrictions t > des fonctions de relaxation ne sont pas toutes les 
fonctions completement monotones, a cause de la condition (|18j) (qui garantit que r est une 
distribution au sens de Schwartz). 

4.2 Relation de conjugaison 

De meme, dans le cas de m parametres observables, les reponses impulsionnelles et les fonctions 
de relaxation sont reliees par : 

(19) = <w 

k 

en notant I la matrice identite m x m. Et par le meme raisonnement, on obtient que les 
fonctions de relaxation sont les distributions de la forme 

(20) r {j (t) = A i:j + B^do + f e~ tx d PlJ (x) 

ou (Aij, (Bij) et (pij) sont des matrices symetriques semi-definies positives, les pij etant des 
mesures telles que : 

r°° 1 

(21) / ——d\ Pij \(x)<+oo. 

J]Q 1 + X 

II resulte de (J2*U|) que les primitives Rij(t) nulles pour t < des r^(t) sont exactement les 
fonctions de la forme (7). Et par consequent sont les reponses impulsionnelles d'un materiau 
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viscoelastique lineaire sans vieillissement que nous appellerons materiau conjugue du materiau 
initial. 

La relation de conjugaison s'exprime sur les reponses impulsionnelles par 
(22) E ^ */*?(*) = ^-fl{t>o}- 

k 

La relation de conjugaison est symetrique et la mise en serie des modeles rheologiques de 
deux materiaux correspond a la mise en parallele des materiaux conjugues et reciproquement. 

Dictionnaires et exemples 

Nous nous limitons ici, pour simplifier, au cas d'un seul parametre observable 



Materiaux 


Reponses impulsionnelles 


processus 


mise en serie 


addition des reponses 
impulsionnelles 


addition de processus 
independants 


mise en parallele 


addition des fonctions 
de relaxation 


addition des Processus 
conjugues independants 


elastique 


f(f\ — -\ r.^T 

- a i {*>0} 


piULcccUo CUIloLcU.lL 

de valeur - 

a 


amortisseur 
(conjugue du 
precedent) 


/(*) = atl {t > 0} 


translation uniforme 
de vitesse a 


materiau de 
Maxwell 


f(t) = (at + b)l {tm 


translation uniforme de 
vitesse a issue du point b 


materiau de 
Kelvin- Voigt 
(conjugue du 
precedent) 


/(t) = I(l-e-f*)l {t > 0} 


processus de Poisson de saut 
d' amplitude | et 
d'intensite - 

a 
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Materiaux 


Reponses impulsionnelles 


processus 


combinaison finie 
d'amortisseurs 
et de ressorts 


/(*) = (Ei " e~ Al *) + at + b)l {tm 


somme de processus de Poissons 
independants de sauts 
d'amplitude A« et d'intensite ai 
et d'une translation de vitesse 
a issue de b 


combinaison finie 
d'amortisseurs et 
de ressorts 
(conjugue du 
precedent) 


f(t) = (EiPi(l-e-> lit )+ct + d)l t >o 

(Ai < Aii < A 2 < H2 < ■ • • si c / 
ou /ii < Ai < H2 < A2 < • • • si a / 0) 


somme de processus de Poissons 
independants de sauts 
d'amplitude \i{ et d'intensite Pi 
et d'un translation de 
vitesse c issue de d 



La relation de conjugaison attire l'attention sur le cas suivant 



Reponses impulsionnelles 


Processus 


/(*) = sfey^o} «e]o,i[ 

conjugue avec le materiau 
d'indice 1 — a 


processus stable unilateral d'ordre a 


seul materiau qui est son propre conjugue 
f(t) = 2/Jl^o} 


processus stable unilateral d'ordre \ 



Autres remarques. 

1 - L'importante propriete du cone des fonctions de Bernstein d'etre stable par composition 
(cf [1]) fait que si fi(t) et /2(f) sont prises dans les tableaux ci-dessus, /i((/2(t)) est encore la 
reponse impulsionnelle d'un materiau viscoelastique lineaire sans vieillissement, ce qui fournit 
une grande variete de fonctions. 

2 - La reponse impulsionnelle d'un materiau viscoelastique a un parametre observable peut 
etre consideree comme une horloge. C'est, en effet, un phenomene declenche a l'instant et 
qui evolue ensuite indefiniment selon une dynamique propre. Ces horloges ont la propriete 
suivante : 



12 



Considerons : 

a) un mobile A anime d'une vitesse rectiligne uniforme, en zero a l'instant 0. 

b) deux curseurs B et C lies respectivement aux reponses impulsionnelles de materiaux (1) 
et (2) et pouvant se deplacer parallelement a A. 

Si on note les positions successives a = 0, a±, • ■ ■ , a n , ■ ■ • de A lorsque B atteint des distances 
regulieres 0, h, 2h, ■ ■ ■ , nh, ■ ■ ■ , alors les positions successives «o = 0, a±, ■ ■ ■ , a n , ■ ■ ■ de A lorsque 
C atteint les a n sont aussi celles qu'il atteindrait lorsqu'un corps viscoelastique D voit son index 
atteindre les distances nh. Cette propriete est due a la stabilite par composition du cone des 
fonctions de Bernstein. 

3 - La theorie de Biot a fait l'objet de nombreuses extensions, d'abord au cas avec vieil- 
lissement, egalement avec des concepts de geometrie differentielle les parametres caches etant 
supposes tensoriels cf [22] et a des cas non lineaires cf [17]. II est vraisemblable qu'une partie 
de ces extensions peut recevoir une interpretation a partir des caracteristiques locales des semi- 
martingales. 
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